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We present a strong coupling expansion that permits to develop analysis of quantum field theory in the infrared limit. Application 
to a quartic massless scalar field gives a massive spectrum and the propagator in this regime. We extend the approach to a pure 
Yang-Mills theory obtaining analogous results. The gluon propagator is compared satisfactorily with lattice results and similarly 
for the spectrum. Comparison with experimental low energy spectrum of QCD supports the view that a resonance is indeed 
a glueball. The gluon propagator we obtained is finally used to formulate a low energy Lagrangian for QCD that reduces to a 
Nambu-Jona-Lasinio model with all the parameters fixed by those of the full theory. 



One of the main difRculties people meets to cope with 
low energy limit of a quantum field theory is the missing of 
a perturbation techniques that could permit to extract an- 
alytical manageable results to compare with experiment. 
In the fall of seventies, Bender's group proposed a possible 
method to solve this impasse |1I2I3I4J . Bender's group ap- 
proach can be understood with a simple quartic oscillator 
as (here and following h = c — 1) 



G(i2-ti)= /[dg(t)]e^/n'^*[^^W^-i^(*)^ 



(1) 



When the limit A — > oo is taken one simply neglects the 
kinematic term assuming it as a higher order correction 
leaving us with 



G(i2-ti)« /[d<z(t)]e-*^n'**^«(*)\ 



(2) 



a highly singular expression that needs to be regularized 
with a proper cut-off. This in turn implies that computa- 
tions produce series generally difficult to resum due to the 
very singular dependence on the cut-off, and the method 
proves to be not effective in the computation of propaga- 
tors and spectra. The reason why Bender's group method 
does not work properly is due to the completely omitted 
dynamics that should be permitted instead. This can be 
accomplished if we do a time rescaling as i ^ ^/\t ~ t 
and we are left with the expression 



G{T2-n)= [dq{T)] 



jVXf:-' drUq{Tf-^,q{rr 



(3) 



that is perfectly regular being just the semiclassical limit 
when A — > oo. In this case the spectrum is computed 
straightforwardly and, accounting for all higher order cor- 
rections, one has the exact result [5]. Starting from the 
Schrodinger equation 
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and finally taking an expansion like ip — "0 + jV' + 
jT^^^^ + O [jj) for the wave function, the above method 
recovers a well-known semiclassical expansion, 
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that is the Wigner-Kirkwood expansion with the spectrum 
given by th usual WKB series [5]. But Wigner-Kirkwood 
expansion is a gradient expansion and this is the key result: 
A strongly coupled quantum system behaves semiclassi- 
cally and its behavior is described by a Wigner-Kirkwood 
gradient expansion. A typical application for this funda- 
mental result could be measurement theory when a mea- 
suring apparatus is strongly coupled to a quantum system. 
But the method is quite general to be applicable to what- 
ever differential equation. 

We would like to extend the above approach to quantum 
field theory. For our aims, we consider a massless quartic 
scalar field theory with a generating functional 



m = / m 



^■,Jd^xU{d<t>f-^<P^+j4 



(7) 



and rescale time as shown above for the quartic oscillator 
to obtain 



that in the limit A ^ oo permits us to recover again the 
semiclassical limit and a gradient expansion [7 . But in this 
case we are able to solve the equation for the propagator 

Giti - ts) + AG(ti - t2f - fi^Siti - t2) (9) 

that is 



(4) G{t,-t2)=e{t,-t2)ii 



sn 



/it 



(10) 



being sn a Jacobi elliptic function and /i an arbitrary inte- 
gration constant. Now, using a small time approximation 
we can take the following definition for the field (p [819] 



m 



dt'G{t~t')j{t') 



(11) 



we can restate the above generating functional for the 
scalar field theory in a Gaussian form 

Z[j] « eH<^Vrf^...^(-v^)5-(..-..).^^^f^.] (^2) 

being 



Zo[j] = 



= pi f d'^^l<i'^^2J(xi)A{xi-X2)j{x2) 



(13) 



with A(a;i -X2) = S^{xi - a;2)[6'(ii - i2)G(ii - ^2) + 0{t2 - 
ti)G{t2 — ti)]. Having obtained the propagator in the 
strong coupling limit we can extract the spectrum of the 
theory in the infrared limit. This can be done using eq. fll 
and the relation for the sn Jacobi function [TOl 



sn(u,i) 



2tt 



oc 



K{i) ^ l + e-(2n+l)7r 



sm 



(2n + 1) 



2K{i) 



being K{i) = / = ^ 

our case for the Feynman propagator 

^1 j-_-[^^„+lg-(„+i)7r 
ri=0 



1.3111028777. This gives in 



AH = ^(2n+l) 



K'^{i) l + e-(2"+i)'^ w2_^2_^^, 



■(15) 



with the harmonic oscillator spectrum 



w„ 



K{i) 



^ ^■ 



(16) 



We see immediately that the propagator is in agreement 
with Kallen-Lehman representation. This is a key obser- 
vation if we want to interpret the poles of the propagator 
as physical states even if the propagator is gauge depen- 
dent. But we also see that this propagator satisfies the 
Callan-Symanzik equation 



,«+«„« + 2.G(0^0 



dpi 



d\ 



(17) 



with a beta function /3 — ~4A and 7 = — 1. This means 
that the theory is trivial in the infrared limit as the cou- 
pling goes to zero as the fourth power of momentum [11 . 

The next step is to apply the above approach to a pure 
Yang-Mills theory in order to analyze the behavior in the 
low energy limit. But we have to cope here with a seri- 
ous difficulty. In the eighties it was proved by Matinyan, 
Savvidy, et al. [12I13I14J that the classical equations of 
motion of Yang-Mills theory are generally chaotic and so, 
completely useless to built a quantum field theory. The 
way out is the following mapping theorem [15j : 

MAPPING THEOREM: An extremum of the action 



S^ 



d^x 



\{d^f - \<^' 



(18) 



is also an extremum of the SU(N) Yang-Mills Lagrangian 
when we properly choose A't with some components being 
zero and all others being equal, and A = Ng^ , being g the 
coupling constant of the Yang-Mills field. 

This theorem holds exactly at classical level but at a 
quantum level can be maintained only in the infrared limit 



as in the ultraviolet case quantum fluctuations spoil the 
mapping producing asymptotic freedom, a property miss- 
ing for the scalar field [16 . Then, using mapping theorem 
we can write down immediately the gluon propagator in 
the Landau gauge as 

/ s 00 2 

DtAv) = ^"(^.^-^)E(2-+l)l^x (19) 



(14) 



p^ 



K^{i) 



(-l)"+ie"("+5)^ 
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1 + e-(2n-l-l)7i- p2_^2_|_j£ 

where we have also introduced Latin indexes for color de- 
gree of freedom and 



"1,1= \ n+ - 



K{z) 



Ng"' 



A. 



(20) 



being A the integration constant. Mapping theorem tells 
us also that the ghost field decouples from the gluon field 
and must behave like a free particle with a propagator 
G{p) = l/(p^ + ie). The running coupling must go to zero 
like the fourth power of momentum being this in agreement 
with lattice computations as shown by Boucaud et al. |17| 
but as also seen from experiments as shown by Prosperi 
et al. |18|19|20] . These latter results on the running cou- 
pling give a strong clue that Yang-Mills theory is trivial as 
happens for the scalar field theory. Lattice computations 
give also results on the propagators [21i22 23^ . These show 
that the ghost is very near the free particle case while the 
gluon propagator reaches a finite value at zero momentum. 
So, we compare our propagator (|19p with lattice computa- 
tions taking the one with the largest volume (27fm)'* [22]. 
For this to work we need to go to Euclidean momenta and 
fix Too in eg. ((20)) . The result is presented in Fig. [l]and is 
really satisfactory. 

The comparison is done in our case with ttiq = 
27^ {Ng'^/2) ■* A « 547 MeV corresponding to a string 

tension ^ = {Ng'^/2)'^A w 456 MeV. This result 
is striking indeed as we observe a mass for the ground 
state very near that of the a resonance and a string ten- 
sion very near that generally obtained from experiments 
410-440±20Mey. So, on the lattice the scenario is clearly 
settled. Numerical solution of Dyson-Schwinger equations 
also confirm it [21] . 

For the spectrum the situation is satisfactory as well. 
Comparing with lattice computations by Teper et al. [25] 
[26] . we have shown in ^7\ that one gets Tab[l] and [2] 
Theoretical numbers compared in the tables are computed 
through the relation ^ = {2n-\- 1)2/^77 that defines a set 
of "golden numbers" that are those measured in lattice 
computations of spectra. 

We cannot say the same with the data by C. Morn- 
ingstar et al. [28] but these authors use anisotropic lat- 
tices. The a state appears on lattice computations of the 
gluon propagator but not on lattice computations of spec- 
tra. It should be said that presently the volumes used in 
the former computations are far larger than for the latter 
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Figure 1. Fit to lattice data of propagator in eq. P^ . 



case. Our view is that an effort should be done in this 
direction in order to get a clear understanding of Yang- 
Mills theory on the lattice by increasing the volumes for 
spectrum computation. 

A needed step is to understand current phenomenology 
of the low energy spectrum of QCD with the above results. 
Indeed, the ground state of a pure Yang-Mills theory can 
be identified with the observed resonance a or /o(600). 
This means that Yang-Mills theory has a higher ground 
state than full QCD as quarks lower it to the pion mass. 
Indeed, /o(600) or a seems to have a large gluonic content 
due to the very small width for 77 decay that can be due 
to TTTT rescattering [29130] and as also seen through xPT 
|31l32j . We have shown that our computed mass is very 
near the one obtained from experiments [33| . If we take 
for the string tension 440 MeV as is customary in lattice 
computations we get ma- = 527MeV and, as seen above, 
lattice computations for the gluon propagator give about 
547 MeV. Mixing with quarks could reduce this value. A 
more conservative choice of 410 MeV can take both a and 
a* , the latter eventually being /o(980), to agree with mea- 
sured data giving nia « 491 MeV and Too-. « 982 MeV. 
This view is shared in recent analysis [29I30J . But while 
these resonances are seen both experimentally and theo- 
retically and, in indirect way through lattice computations 
for the propagator, a direct evidence on lattice computa- 
tions for spectra appears difficult to achieve presently [34| . 

As pointed out in [3S], once gluon propagator is known 
one can obtain a low energy model directly from QCD. We 
derived this model in [36] . So, taking QCD action as 



Excitation 


Lattice 


Theoretical 


Error 


a 


- 


1.198140235 


- 


0++ 


3.55(7) 


3.594420705 


1% 


0++* 


5.69(10) 


5.990701175 


5% 



Table 1 
Comparison for the 0++ glueball 



S = Jd^x Y,q(n-d-gYl-A''-m,)q (21) 



1g';^^g^^'' - ^{d ■ Af 



being G°^ = 9^yl° - d^A''^+ gf'''"'AlAl, we use the small 
time approximation seen for the quartic scalar and Yang- 
Mills fields [819] with the above computed propagator writ- 
ing 



A^^ix) « g I d'yDf,{x - y)J2 i{y)^Yq{y) (22) 



Excitation 


Lattice 


Theoretical 


Error 


cr* 


- 


2.396280470 


- 


2++ 


4.78(9) 


4.792560940 


0.2% 


2++* 


- 


7.188841410 


- 



Table 2 
Comparison for the 2++ glueball 



givmg 



^q{x){i-i ■ d ~ ■mq)q{x) 



(23) 



L q 



-g'^^q{x)—Yq{x) X 



9:9' 



d'yDf^{x-y)q'{v)-Yq'{y) 



In the infrared limit p <C mo, being toq = 
{Ng'^/2)iTr/2K{i)A = m^ « A91MeV that can be taken 



as the NJL cut-off Anjl, one gets 



d^x 



^q{x){ij ■ d ~ mg)q{x) 



(24) 



:,g^^q{x)—-f''q{x) X 



d'^yGix - y)q'{y)—i^,q'{y) 

being G{x — y) ~ 3.761402959 • ■^S'^{x — y) (Fermi approx- 
imation). This gives in the end 



S 



d^3 



^q{x){i-fd' 



M^) 



(25) 



^Gnjl Y^ q{x)—l^q{x)q'{x)—-i^q'{x) 



being Gnjl ~ 3.761402959 • ^ that links QCD theory 
parameters and the couphng of the NJL model. The point 
opened up by the introduction of quarks is how to recover 
gluonic contributions to QCD spectrum. This possibility 
can be achieved through bosonization techniques [37] . 

So, a consistent formulation of quantum field theory in 
the infrared exists. We have seen that a satisfactory agree- 
ment can be obtained in this way with lattice computa- 
tions both for the propagators and the spectrum for Yang- 
Mills theory. But while for the propagators the agreement 
is really good, for the spectrum some more effort on lat- 
tice computations seems needed to obtain a satisfactory 
comprehension of all the matter. From the experimental 
side the situation is quite good granting an understand- 
ing of where glueballs should lie in the observed spectrum. 
Agreement is seen between different theoretical approaches 
in this case while some more work is needed to clarify 
the situation. It should be said that if further elements 
will confirm this scenario, surely a lot of unexpected views 
about strong interactions will force new understanding for 
all quantum field theory. 

REFERENCES 



1. 



2. 



3. 



4. 
5. 

6. 

7. 



C. M. Bender, F. Cooper, G. S. Guralnik, and D. H. 
Sharp, Phys. Rev. D 19, (1979) 1865 . 
C. M. Bender, F. Cooper, G. S. Guralnik, D. H. Sharp, 
R. Roskies, and M. L. Silverstein, Phys. Rev. D 20, 

(1979) 1374 . 

C. M. Bender, F. Cooper, G. S. Guralnik, H. Moreno, 
R. Roskies, and D. H. Sharp, Phys. Rev. Lett. 45, 

(1980) 501 . 

C. M. Bender, F. Cooper, R. Kenway, and L. M. Sim- 
mons, Phys. Rev. D 24, (1981) 2693 . 
C. M. Bender, K. Olaussen and P. S. Wang, Phys. Rev. 
D 16, (1977) 1740 . 

M. Frasca, Proc. R. Soc. A (2007) 463, 2195. 
M. Frasca, Phys. Rev. D 73, (2006) 027701 ; Erratum- 
ibid., (2006) 049902. 



8. M. Frasca, Mod. Phys. Lett. A 22, (2007) 1293. 

9. M. Frasca, Int. J. Mod. Phys. A 23 (2008) 299. 

10. I. S. Gradshteyn, I. M. Ryzhik, Tahlc of Integrals, Se- 
ries, and Products, (Academic Press, 2000). 

11. M. Frasca, ■arXiv:0802. 1183" [hep-th]. 

12. S. G. Matinyan, G. K. Savvidy, N. G. Ter-Arutunian 
Savvidy, Sov. Phys. JETP 53, 421 (1981). 

13. G. K. Savvidy, Phys. Lett. B 130, 303 (1983). 

14. G. K. Sav vidy, Nucl. Phys. B 246, 302 (1984). 

15. M. Frasca. [arXiv:0709.2042l [hep-th] . 

16. M. E. Peskin, D. V. Schroeder, An Introduction to 
Quantum Field Theory, (Perseus, Reading, 1995), 
Ch.l2. 

17. Ph. Boucaud, F. De Soto, A. Le Yaouanc, J. P. Leroy, J. 
Michcli, H. Moutarde, O. Pene, J. Rodrguez-Quintero, 
JHEP 0304 (2003) 005. 

18. G. M. Prosperi, M. Raciti, C. Simolo, Prog. Part. Nucl. 
Phys. 58, 387 (2007). 

19. M. Baldicchi, A. V. Nesterenko, G. M. Prosperi, D. 
V. Shirkov, C. Simolo, Phys. Rev. Lett. 99, 242001 
(2007). 

20. M. Baldicchi, A. V. Nesterenko, G. M. Prosperi, C. 
Simolo, arXiv:0705.1695' [hep-ph]. 

21. I.L. Bogolubsky, V.G. Bornyakov, G. Burgio, E.M. II- 
genfritz, V.K. Mitrjushkin, M. Miiller-Preussker, P. 
Schemel, PoS (LATTICE 2007) 318. 

22. A. Cucchieri, T. Mendcs, PoS (LATTICE 2007) 297. 

23. A. Sternbeck, L. von Smekal, D. B. Leinweber, A. G. 
WiUiams, PoS (LATTICE 2007) 340. 

24. A. C. Aguilar, A. A. Natale, JHEP 0408 (2004) 057. 

25. B. Lucini, M. Teper, Phys. Rev. D 64, 105019 (2001). 

26. B. Lucini, M. Teper, U. Wenger, JHEP 06, 012 (2004). 

27. M. Frasca, arXiv:0704.3260 [hep-th]. 

28. Y. Chen, A. Alexandru, S. J. Dong, T. Draper, I. Hor- 
vath, F. X. Lee, K. F. Liu, N. Mathur, C. Morningstar, 
M. Peardon, S. Tamhankar, B. L. Young, J. B. Zhang, 
Phys. Rev. D 73, 014516 (2006). 

29. G. Mennessier, S. Narison, W. Ochs, 0804.4452 [hep- 
ph]. 

30. G. Mennessier, P. Minkowski, S. Narison, W. Ochs, 
0707.4511 [hep-ph]. 

31. J. R. Pelaez, G. Rios, Phys. Rev. Lett. 97 (2006) 
242002. 

32. C. Hanhart, J. R. Pelaez, G.Rios, 0801.2871 [hep-ph]. 

33. D.V. Bugg, 0804.3450 [hep-ph]. 

34. C. McNeile, 0710.2470 [hep-lat]. 

35. T. Goldman, R. W. Haymaker, Phys. Rev. D 24, 724 
(1981). 

36. M. Frasca EXh.: 0803.0321 [hep-th]. 

37. D. Ebert, [hep^/9710511| 



